Novel two-loop SUSY effects on CP asymmetry in B — > <pK t 
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Inspired by the exotic measurements on the CP asymmetry in B — » 4>K S , we study a new diagram 
in supersymmetric models which can make the difference sin 24>\ s 1 {J/<HK S ) - sin 2(j>\ 1! (</>K 3 ) to be 
20 — 50% after satisfying the constraint from 6 — > S7. We also find that the direct CP asymmetry 
of b — > S7 could be ~ 10% and testable at B factories. 
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While enjoying the large CP asymmetry (CPA) in the 
decay of B — ► J/^K S observed by Belle [lj and Babar 
at the precision level, the recent data on B — ► tt + tt~ |3j 
and B -► ^ 0,0 have stimulated theorists to think 
more about other possible CP violating phases, beside 
the Kobayashi-Maskawa (KM) [(| phase in the standard 
model (SM). 

It is known that with the Wolfenstein parametriza- 
tion |7j, the tree and penguin diagrams have the same 
CP phase for the inclusive processes of b — > sec and 
b — > sss. Thus, the time-dependent CPA, proportional 
to T(B -> f C p)~ r(B -> /cp) with / CP being the final 
state and having a definite CP property, arises from the 
B — B oscillation dictated by box diagrams, in which the 
source of the CP phase is from V td = | V^^ | e~ £ ^ >1 . For the 
channel of fcp = J/^K S , the CPA is related to sin20i 
and the mixing-induced CP violation. If there is only 
the KM phase involved in the low-energy, the pure pen- 
guin process of B — > 4>K S has approximately the same 
value of sin 20i as that in the decay of B — > J/^>K S , i.e., 
AS* 01 = sin20i(J/* J ftT s ) - sin 2<j) 1 (<t>K s ) ~ 0. 

It is usually believed that new physics could go into 
low-energy phenomena through loop diagrams, in which 
new particles appearing in the loops are integrated out 
and the remaining effective couplings are as functions of 
their masses and couplings to the conventional particles. 
Since the transition of b — > sss is a pure quantum loop 
effect, one can recognize immediately that B — ► 4>K S is 
a good candidate to probe new physics. Furthermore, 
although the tree-level contributions in b — > sec are over 
a factor of 5 larger than those of penguin diagrams in 
the SM 9], the penguin-type diagrams induced by new 
physics could be enhanced, which will clearly affect the 
decay of B — > J/^K S , especially on its direct CPA. 

To understand the Belle's result of the 3.5er difference 
on sin 20i between J/^K S and cf>K s modes 0, vari- 
ous theoretical models such as those with supersymme- 
try (SUSY) [H EE 13 and left-right symmetry 
have been investigated. In addition, the authors of Refs. 
[ToL ITU have tried to solve the problem of unexpected 
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large branching ratios (BRs) in B — > rj'K decays. How- 
ever, we would like to address some problems on these 
attempts as follows: 

(a) Direct CP violation on B — > J/^K S : 

We emphasize that Belle and Babar not only measure 
an accurate mixing-induced CPA, but also indicate no 
direct CPA in B — > J /^K s , up to the percentage level. 
Those new SM-like effective four-fermion interactions for 
b — > sss will inevitably contribute to b — > sec. It is also 
known that there exist large strong phases in the produc- 
tion of charmed mesons (including charmonium states) 
[H Therefore, to enhance the BRs of B -> rj'K s 
with large CP violating effects will make the direct CPA 
in B — > J/^K S to be over the current experimental lim- 
its. 

(b) BRs ofB^nK and B -> rf'^K*: 

We note that the problems for the production of 77' 
in B decays depend on not only B — * ij'K, but also 
B — * i]K and B — > jf'^K*. From the data at Babar, 
we have that BR(B -> r]K a ) = (2.9 ± 1.0 ± 0.2)1(T 6 
[13, BR(B -> ry.fi:* ) = (18.6 ± 2.3 ± 1.2)1CT 6 , and 
BR(B -> 77'^*°) < 7.6 x 1CT 6 0. By using the pertur- 
bative QCD approach ^^|, we find that the estimating 
BRs of B — > ryif and t]'K*° are over the current exper- 
imental values, whereas it is lower for £? — > r/K* . 

It is clear that to resolve the problems we need more 
knowledge on 77^ mesons as well as their relevant physics. 
On the other hand, we may bypass these problems by 
concentrating on new physics effects which are insensi- 
tive to hadronic uncertainties. In this paper, we will 
introduce a two-loop diagram illustrated in Fig. ^ in 
the framework of SUSY models, resulting from dipole 
operators. In contrast with other mechanisms, such as 
those discussed in Refs. 0,^3 m which the relevant off- 
diagonal terms of squark-mass matrices directly involve 
flavor changing neutral current that couples to gluino, 
our two-loop effect shows how to generate the flavor 
changing processes naturally in the SUSY models. We 
will illustrate that the diagram not only contributes a 
sizable value for the difference of sin20i between J/^K S 
and 4>K S channels, but also satisfies the experimental 
constraints such as those from the b — > s"f decay and 
the neutron electric dipole moment (NEDM). Since the 
diagram involves the couplings of the charged Higgs to 
squarks, we first discuss the relevant couplings in SUSY 
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FIG. 1: Two-loop diagram by the ba — §l flavor changing 
effect and the chromodipole operator. 

models, given by pcj 

C Hff = -{2V2G F y! 2 (y tb A b L tlb R 

+V t * s A R t* R ~s L )H+ + h.c, (1) 

where A b L = m b (Al tan/3 — /x) and A R = mt(A t cot/3 — 
/x*). Here, the definition of the angle (3 is followed by 
tan/3 = v u /vd with w M and being the vacuum ex- 
pectation values (VEVs) of Higgs fields $ u and <E> d re- 
sponsible for the masses of upper and down type quarks, 
respectively, and \i is the mixing effects of $"> d . For 
a large tan/3 case, A\ and A R can be simplified as 
A\ w mf,Af,tan/3 and w — m t /i*. Note that we have 
neglected the contribution of sr because the correspond- 
ing coupling is associated with the strange-quark mass. 
Moreover, in order to suppress one-loop contributions, we 
assume that the flavor mixing effects on the down-type 
squark mass matrix are small. 

We remark that both flavor changing and chirality flip- 
ping are involved in Fig. ^ in which the charged Higgs 
is used to change the flavor and the mixing of fx and t R 
to govern the chirality flipping, representing by the cross 
in the figure. Explicitly as usual, the mixing terms are 
described by |2l| 

= (M^-Zicot^mu, (2) 

where (M^)i^ represent the trilinear soft breaking ef- 
fects. For simplicity we have adopted the so-called super- 
CKM basis, where quarks are in the mass eigenstates so 
that mu is the diagonal upper-type quark mass matrix 
|2l||. To overcome the NEDM constraint, it has been 
proposed [2^| to use hermitian Yukawa and A matrices. 
The construction of a hermitian Yukawa matrix can be 
implemented based on some symmetries, such as the hor- 
izontal SU(3)h H3 and left-right |24| symmetries. As a 
result, the CP phases of O(l) can exist naturally even 
with the NEDM contributions. Moreover, it implies that 
the CP asymmetries in hyperon decays could reach the 
value of O(10 -4 ) [2{|, which is testable in the experi- 
ment E871 at Fermilab [2i|. However, in the class of 



models proposed in Ref. [22j, the /x parameter is real 
which is not favored in our following discussions. To 
avoid this shortcoming, we address the NEDM constraint 
by imposing the Yukawa and A matrices to be hermi- 
tian and the squark mass of the first generation to be 
O(10) TeV. Hence, the /x parameter is regarded as a com- 
plex value in our approach. Due to the hermitian prop- 
erty, a special relation is obtained as \&ki) LR = \&ki) RL 
with (8%) LS = (M 2 0kl ) LR /m 2 = (V u ^A u U u V u ) kl /m 2 , 
where A u ^ — A u , V u is the mixing matrix for diago- 
nalizing the mass matrix of upper-type quarks and m 
is the average squark mass in the super-KM basis. In 
general, the trilinear SUSY soft breaking A® terms are 
not diagonal matrices. However, due to the relation 
of A% = (Y Q A Q )ij with Y Q (A Q ) being Yukawa (A- 
parameter) matrices and the small effect of renormal- 
ization group, dominant effects of A® are still from the 
diagonal elements j2^ if we take A® to be universal and 
diagonal at the grand unified scale. We use Aq — A® to 
simplify our estimations. Therefore, the contribution in 
Fig. ^is proportional to m b m t ^* A b tan (3(S 33 )l R - Since 
A u is hermitian, A t (5 f 33 )LR can be regarded as real val- 
ues. Hence, in our mechanism, the CP violating source 
is focused on the complex /x term. We note that by 
adopting a large tan/3, the /x-dependent effect is from 
the vertex of the charged Higgs coupling to squarks. For 
convenience, we write the relationship between weak and 
physical eigenstates for the mixing of ti and t R as 

f t L \ _ / cost9 t sin 6 t \ f h\ /o\ 
\t R J \-sm9 t cos6 t )\t 2 )' [ ' 

To study Fig. ^ we start with the effective interactions 
for quark-gluino-squark, given by |27| 

£gq q = -V2g s {sP R g a T a S L - bP L ~g a T a b R ) + h.c, (4) 

where the flavor mixings for squarks have been neglected. 

It is interesting to note that if we use the photon in- 
stead of the gluon and include the emission of the photon 
at the charged Higgs, we find that the same mechanism 
could also contribute to b — » 57. Therefore, sizable val- 
ues for both AS c j >1 and the rate CPA in B — > X s -f can 
definitely provide a hint for new physics. The effective 
operators for b — > sj(g) are given by 

C = ^V t * s V tb (C 77 (fi)0 7l + C Sg ([i)0 Sg ), (5) 

where 77 = m b e/(8n 2 )sa^F^(l + 75)6, O gg = 
m b9s /(8ir 2 )sa^T a G a ^{l + 75)6, 
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C Sa = C 7l /{2N C PJ) and Pj = C F (Qu - 1) with C F = 
4/3 and Qu — 2/3 being the color factor and the charge 
of the upper-type squark, respectively. Clearly, we obtain 
the unique property that the effects of electric and mag- 
netic dipole moments are directly related to those of chro- 
moelectric and chromomagnetic dipole moments, respec- 
tively. Before we proceed further, we have to examine 
whether the two-loop effects are of interest. Explicitly, we 
would like to check whether the value of CV 7 is larger or 
smaller than experimental constraint 0.3 < \C?~\ < 0.34 
|29| . For an illustration, we set the values of parameters, 
by satisfying the constraints from the NEDM [3(J , as fol- 
lows: tan/3 ~ m t /m bl sin6> f cos6» t ~ 0.2, V t *V tb /V t *V tb ~ 
O(l), A b /rrig ~ /i/mg ~ 4, Arg{^) = tt/2, m# ~ 150 
GeV, m g ~ 1 TeV, m t ~ ~ 200 GeV and m b = m s ~ 500 
GeV, and we have |CV 7 | ~ 0.8. If we take tan/3 ~ 50, 
sin 6^ cos # t ~ 0.35, A b /m g ~ /V m s ~ 0(1) and the re- 
mains to be the same as the above choices, we obtain 
|CV 7 | ~ 0.14. Furthermore, by using b-quark and sbot- 
tom instead of s-quark and its squark, the similar two- 
loop diagram could contribute to the EDM of s-quark. 
It is known that the current limit of the s-quark chromo 
EDM is |edf \ expt < 5.8 x 10~ 25 e cm 3! . We now exam- 
ine the contribution to in our mechanism. By using 
Eq. (J5J and assuming m^ = m^ and A s = A b , we obtain 



V2G F V 2 ^\ImC 8g \ = 



GfV 2 



V2C F 8tt 



-Im\G 
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(6) 



Numerically, we get |ed~ | ~ 2.5 x 10 _25 |CV 7 | e cm, which 
is below |ed^| e2 , p t if |CV 7 | < 1. Clearly, in our mechanism 
it is inevitable to utilize the large tan [3, A b jm- g and /x/m§ 
scheme and, therefore, the most strict constraint is the 
BR of B -> X sl . 

In order to discuss the mixing-induced CP problem in 
B — > (j)K s , we write the relevant definition of the time- 
dependent CPA as 



A C p = 



BR(B -+<j>K a )- BR(B -> <f>K a ) 



BR(B -> (j)K s ) ■ 
Crt,K s cos Amst - 

|A|2 " 1 A * 

cos /Xmst 



BR(B -> 
S^k, sin Amst, 
21mA 
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sinAmsi, (7) 



where A = fT ^ K -)A(B -> cpK s )/A(B -> 0Jf a ) and 

— > /cp) is the decay amplitude. Since the dipole 
operators contributing to the nonleptonic decays belong 



to next-to-leading order in a s , we can safely neglect the 
contributions to the decay amplitude of B — > J/^>K S . 
For displaying the other SUSY effects on the B — B mix- 
ing, we use <^j^ instead of fa. Hence, 4>\^ is still deter- 
mined by B — ► J/^K Sl exclusively. For estimating the 
hadronic matrix element of _B — + </>if, we use the naive 
factorization, given by 



(cf>K\0 Sg \B,p B ) 



2a s ml 
9ir 



■Pb,{8) 



where F BK (Q) is the transition form factor of B — > K 
at Q 2 = 0, q 2 is the squared momentum of the virtual 
gluon, e, f$ and m^ correspond to the polarization vec- 
tor, decay constant and the mass of </>, respectively. The 
dominant contribution of factorization assumption is con- 
firmed by the PQCD approach [32| in which q 2 is related 
to the momentum fractions of quarks and convolutes with 
wave functions. We note that although 07 7 can also con- 
tribute to the decay of B — > faK s , since the coupling is 
electromagnetic interaction and much smaller than that 
of strong interaction, we neglect its contribution. Ac- 
cordingly, the decay amplitude for B — > 4>Kq is written 
as 



A{B^4>K ) = ~qV t * s V tb ]>> 



9tt q 2 



xf 4> m 4> F BK (Q)e* -p B , 



(9) 



where dj, defined in Ref. [33], stand for the effective 
Wilson coefficients in the SM, included from electromag- 
netic penguin diagrams. The value of 5Zi=3 a * ^ s es ti- 
mated to be —0.045. The parameter A in Eq. J7J for the 
CPA can be simplified as A = e -^<l>f(J/^K,) e -i2<t, New = 
e -te4>f(<t>K a ) with 
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(10) 
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To display the unique character of the two-loop di- 
agram, we adopt the value of C7 7 such that C7 7 = 
—Cf^ 1 ± i\ImC^ | and the experimental value C*l* = 
+ C*7 7 = ±z|ImC7 7 | instead of scanning the whole 
parameter space. By using C^ 1 = —0.30 and the iden- 
tity Cs g = —3Cj 1 /8, the CP violating phase from the de- 
cay amplitude is tan favew = T(0.18±0.0ljlo e)j in which 
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the first error is from |Cyr| = 0.32 ± 0.02 and the sec- 
ond theoretical error arises from the uncertainty in q 2 = 
(3/8 ± l/8)m|. Since S^ Ks = sin2<pf(<f>K s ), by taking 
sin2</>f T (J/ , I>iC j ) sa 0.74 measured by Belle and Babar, 



we obtain S 4>Ks = 0.46±0.01I^(0.93±0.01+^) where 
the sign of (f>New is chosen to be negative (positive). In- 
terestingly, the former value is close to the central value 
of the Babar's result p. Furthermore, we can straight- 
forwardly calculate the difference of the CPAs to be 



A,eff 



sin2(j)l ff {J/^K s ) - sin2(t)l ff (<j>K s ) 
0.28±0.0lt°;?£ (-), 



-(0.20 ±0.011°;°*) (+) 



(11) 



We now consider the two-loop effects for the CPA in b — > 
sj. According to the formalism shown in Ref. j^S], the 
rate CPA for b — > sj is given by 



Acp (b -> sj) 



1 



100 |Cfff| 

+9.52ImC^C|f 



{l.lImC 2 Cff 



0.16Ima 2 C 8 e f }, 



where C 2 » 1.11 and = C§ g M + C$ g . With 

the same CV 7 used above, we get Acp(b — > 57) » 
±(10.5±0.6)% for negative and positive signs in ImC7 7 , 
respectively. Comparing to the recent Babar's limit of 
-0.06 < A C p(b -> S7) < +0.11 jlHl, we find that only 



the result with negative sign in Eq. (|llfl is reliable, which 
could be used to resolve the sign ambiguity in ImC7 7 . Fi- 
nally, we remark that although our upper value on the 
CP asymmetry of 6 — > S7 is a little bit over the Babar up- 
per bound, the problem can be removed by relaxing the 
required condition CV 7 = —C^ ± illmC^! introduced 
for our simplified analysis. 

In summary, we have studied the novel two-loop SUSY 
effects on the CPAs of B — > (j)K 8 and b —* sj. We have 
found that with large values of tan/? and A^{ji)/mg, the 
difference of sin2^ ff between J/^K S and <pK s can have 
a deviation of 20 — 50%. The main theoretical error is 
due to the uncertainty in q 2 . We have also shown that 
the two-loop effect can give the CPA in b — > sj around 
+10%. It is clear that, since the two-loop contributions 
to the CPAs in both decay modes can be the dominant 
ones in the SUSY models, experimental measurements 
at B factories on these CPAs can determine the sizes of 
these novel contributions. 

Note added: Our two-loop SUSY mechanism has been 
applied to the decay of B s — > n + ^~ [3fij . 



Acknowledgments 

This work is supported in part by the National Science 
Council of R.O.C. under Grant No. NSC-91-2112-M-001- 
053, No. NSC-92-2112-M-006-026andNo. NSC-92-2112- 
M-007-025. 



Belle Collaboration, K. Abe et al, Phys. Rev. D66, 
071102 (2002). [21] 
Babar Collaboration, B. Aubert et al, Phys. Rev. [22] 

Lett. 89, 201802 (2002). [23] 

Belle Collaboration, K. Abe et al, hep-ex/0401029 [24] 
Belle Collaboration, K. Abe et al, Phys. Rev. Lett. 91, [25] 

261602 (2003). [26] 

Babar Collaboration, B. Aubert et al, hep-ex/0403026 [27] 
M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, [28] 
652 (1973). 

L. Wolfenstein, Phys. Rev. Lett. 51, 1945 (1983). [29] 
Y. Grossman and M.P. Worah, Phys. Lett. B395, 241 
(1997); D. London and A. Soni, Phys. Lett. B407, 61 
(1997). 

Y.H. Chen et al, Phys. Rev. D60, 094014 (1999). 

S. Khalil and E. Kou, Phys. Rev. Lett. 91, 241602 (2003). [30] 

B. Dutta et al, Phys. Rev. Lett. 90, 011801 (2003). [31] 

C. K. Chua et al, Phys. Rev. Lett. 92, 201803 (2004). 

A. Datta, Phys. Rev. D66, 071702 (2002); J.F. Cheng et [32] 
al, Phys. Lett. B585, 287 (2004). 

M. Raidal, Phys. Rev. Lett. 89, 231803 (2002). [33] 

C.H. Chen, Phys. Rev. D67, 094011 (2003). [34] 

Y.Y. Keum et al, Phys. Rev. D 69, 094018 (2004). 

Babar Collaboration, Phys. Rev. Lett. 92, 061801 (2004). [35] 

Babar Collaboration, B. Aubert et al, hep-ex/0403025 

E. Kou and A.I. Sanda, Phys. Lett. B525, 240 (2002). [36] 

J.F. Gunion et al, The Higgs Hunter's Guide (Addison- 



Wesley, Reading, M.A., 1990). 

A.J. Buras et al, Nucl. Phys. B520, 3 (1998). 

S. Abel et al, Phys. Lett. B504, 241 (2001). 

A. Masiero and T. Yanagida, hep-ph/9812225] 

K.S. Babu et al, Phys. Rev. D61, 091701 (2000). 

C.H. Chen, Phys. Lett. B521, 315 (2001). 

C. White et al, Nucl. Phys. Proc. Suppl. 71, 451 (1999). 

H.E. Haber and G.L. Kane, Phys. Rept. 117, 75 (1985). 

S. Bertolinia and F. Vissani, Phys. Lett. B324, 164 

(1994). 

ALEPH Collaboration, R. Barate et al, Phys. Lett. 

B429, 169 (1998); Belle Collaboration, K. Abe et al, 

Phys. Lett. B 511, 151 (2001); CLEO Collaboration, S. 

Chen et al, Phys. Rev. Lett. 87, 251807 (2001); BaBar 

Collaboration, B. Aubert et al, |hep-ex/0207076| 

A. Pilaftsis, Nucl. Phys. B644, 263 (2002). 

J. Hisano and Y. Shimizu, Phys. Rev. D70, 093001 

(2004); G.L. Kane et al, hep-ph/0407351 

S. Mishima and A.I. Sanda, Phys. Rev. D69, 054005 

(2004). 

C.H. Chen et al, Phys. Rev. D64, 112002 (2001). 
A.L. Kagan and M. Neubert, Phys. Rev. D58, 094012 
(1998); K. Kiers et al, Phys. Rev. D62, 116004 (2000). 
Babar Collaboration, B. Aubert et al, Phys. Rev. Lett. 
93, 021804 (2004). 

S. Baek, Phys. Lett. B595, 461 (2004). 



